Correlation effects in ionic crystals: 

*o ! I. The cohesive energy of MgO 

o\ '■ 

1 : Klaus Doll, Michael Dolg, Peter Fulde 

a ■ 

^ ■ Max-Planck-Institut fiir Physik komplexer Systeme 

. D-01187 Dresden, Germany 

Hermann Stoll 

On 

O Institut fiir Theoretische Chemie 

| Universitat Stuttgart 

O ' D-70550 Stuttgart, Germany 

in : 

February 1, 2008 

a 

Abstract 

U ; 

Ch , High-level quantum-chemical calculations, using the coupled-cluster 

approach and extended one-particle basis sets, have been performed 
for (Mg 2+ ) n (0 2 ~) m clusters embedded in a Madelung potential. The 
results of these calculations are used for setting up an incremental ex- 
, pansion for the correlation energy of bulk MgO. This way, ~ 96% of 

' the experimental cohesive energy of the MgO crystal is recovered. It is 

shown that only ~ 60% of the correlation contribution to the cohesive 
energy is of intra-ionic origin, the remaining part being caused by van 
der Waals-like inter- ionic excitations. 
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1 Introduction 



While the density-functional (DFT) method, with its long tradition in solid- 
state physics, is getting wide acceptance in the field of quantum chemistry 
nowadays, there are also attempts to the reverse, i.e., making use of the tra- 
ditional quantum-chemical Hartree-Fock (HF) and configuration-interaction 
(CI) methods not only for molecular but also for solid-state applications. 
The Torino group of Pisani and co-workers, e.g., devised an ab-initio HF 
scheme for solids [1], which has successfully been applied to a broad range 
of (mostly) covalently bonded and ionic solids, within the past five years. 
A main asset of the HF scheme is the availability of a well-defined wave- 
function, which may not only be used for extracting properties but also as a 
starting-point for systematically including electron-correlation effects. Such 
effects, which are only implicitly accounted for in density- functional meth- 
ods, often have a strong influence on physical observables, in molecules as 
well as in solids. Several suggestions have been made how to explicitly in- 
clude electron correlation in solids, among them the Quantum Monte-Carlo 
(QMC) approach [2], the Local Ansatz (LA) [3], and the method of local 
increments [4] (which may be considered as a variant of the LA); in QMC 
the HF wavefunction is globally corrected for electron-correlation effects by 
multiplying it with a factor containing inter-electronic coordinates (Jastrow 
factor); the latter two methods rely on applying selected local excitation 
operators to the HF wavefunction and thus have a rather close connection 
to traditional quantum-chemical post-HF methods. 

The number of test examples is still rather limited with all three solid- 
state correlation schemes, and is mainly restricted to semiconductors so far. 
For ionic insulators where quantum-chemical methods would seem to be 
most suitable and easily advocated, much work indeed has been devoted to 
correlation effects on band structures (cf. e.g. [5a,b]), but only few stud- 
ies refer to cohesive energies (cf. e.g. [5c]), and only a single application of 
the post-HF schemes mentioned in the last paragraph exists to our knowl- 
edge (NiO with QMC [5d]). This does not mean that such applications to 
ground-state correlation effects are without challenge. For MgO, the system 
to be dealt with in this paper, HF calculations [6] yield a lattice constant 
which is in agreement with experiment to ~ 0.01 A, but the correlation con- 
tribution to the cohesive energy is significant (~ 3 eV, nearly half as large 
as the HF value). The local-density approximation (LDA) of DFT does 
not a good job here either: an overbinding results which is more than two 
times as large as the correlation contribution to the lattice energy [7], and 
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invocation of gradient-corrected functionals is indispensable for obtaining 
reasonable results [6b], cf. Sect. 3.5 below. The situation would not seem 
too complicated, nevertheless, if the effect could be explained just by adding 
correlation contributions of individual ions; in fact, such a suggestion has 
repeatedly been made in literature, cf. Refs. [6a, 17c]. However, the O 2- ion, 
one of the building blocks of the MgO crystal, is not stable as a free entity, 
and an accurate determination of the correlation energy of this highly fluc- 
tuating and easily polarizable ion in its crystal surroundings is not expected 
to be an easy task; the more so, since already for the determination of the 
electron affinity of the free O atom high-level quantum-chemical correlation 
methods are required [8]. Moreover, as we shall show below, intra-ionic in- 
teractions contribute only with ~ 60% to the total correlation effect on the 
bulk cohesive energy, and van der Waals-like inter-ionic excitations play an 
important role here. 

This paper is the first in a series devoted to the application of the method 
of local increments to ionic solids; it shows, at the example of MgO, how 
to set up an incremental expansion of the bulk correlation energy using in- 
formation from quantum-chemical calculations on finite clusters only. The 
organization of the paper is as follows: in Section 2, computational de- 
tails are given for the applied quantum-chemical methods, and test calcula- 
tions are performed for the first two ionization potentials of the Mg atom, 
for the electron affinity of the O atom, and for spectroscopic properties of 
the MgO molecule. The method of local increments is briefly described in 
Section 3; correlation-energy increments are evaluated from calculations on 
(Mg 2+ ) ra (0 2_ ) m clusters, and the incremental expansion for the total cor- 
relation energy of bulk MgO is discussed. Conclusions follow in Section 
4. 

2 Test calculations 

Our first test concerns the electron affinity, EA, of the oxygen atom. We 
performed separate calculations for the ground states of O and 0~, respec- 
tively, using various single- and multi-reference quantum-chemical config- 
uration interaction methods for treating many-body correlation effects [9]. 
The one-particle basis sets employed have been taken from the series of 
correlation-consistent (augmented) polarized Gaussian basis sets of Dun- 
ning and co-workers [10,8a]. The ab-initio program package MOLPRO [11] 
has been used in these and all of the following calculations of the present 
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paper. The results for EA(O) are collected in Table 1. It is seen that the 
best single- and multi-reference methods (coupled-cluster with single and 
double excitations and perturbative inclusion of triples (CCSD(T)) and 2s - 
2p active space; multi-reference averaged coupled pair functional with single 
and double excitations (MR-ACPF) and 2s - 3p active space) yield quite 
similar results (1.40 eV), which differ from the experimental value by only 
0.06 eV. Effects of the one-particle basis set are significant, even at the stage 
of including g functions, and probably are responsible for the major part of 
the remaining deviation from experiment. In the (Mg 2+ ) n (0 2 ~) m cluster 
calculations to be described in Sect. 3, we could only afford the valence 
triple-zeta [5s4p3d2f] basis set, at the single-reference level; the concomi- 
tant differential errors can be estimated to about 0.1 eV per O atom and 
(added) electron. 

The next test deals with the ionization potentials, IP, of the magnesium 
atom (Mg — > Mg + , Mg + — > Mg 2+ ). A difficulty is encountered here, since 
not only the correlation energy of the valence (3s 2 ) electron pair has to 
be accounted for but also core-valence correlation effects are non-negligible: 
the latter contribute with 0.3 eV to the Mg + — > Mg 2+ IP, e.g.. For ac- 
curately describing these effects explicitly, a high computational effort is 
needed in quantum-chemical ab-initio calculations. In (all-electron) calcu- 
lations with a basis set of medium quality ((12s9pld)/[5sApld\ [12]) errors 
of 0.12 and 0.21 eV remain for the two IPs, and even a very large uncon- 
tracted (20sl5p6d3f) basis [13] still yields deviations from experiment of 
0.03 and 0.04 eV in CCSD(T) calculations. Without loss of accuracy, how- 
ever, the computational effort can be effectively reduced [14,15] by simulat- 
ing the Mg 2+ core by a pseudopotential (PP) which describes core-valence 
interaction at the HF level, in conjunction with a core-polarization poten- 
tial (CPP) which accounts for core-valence correlation effects. Using these 
methods, very good agreement with experiment is obtained, at the corre- 
lated level, cf. the results of Table 2. (Only ACPF data are given, in the 
Table, since for a two-electron system all the correlation methods of Table 
1 (CI, ACPF, CCSD) coincide.) For the cluster calculations of Sect. 3, we 
adopt the PP+CPP description, together with the energy-optimized (4s4p) 
valence basis set; the concomitant differential errors can be estimated to 
about 0.02 eV per Mg atom and (removed) electron. 

The last test of our quantum-chemical arsenal of methods was performed 
for the MgO molecule. The basis sets applied here are the same as those 
used in the next section for (Mg 2+ ) n (0 2_ ) m clusters. The results for bond 
length, dissociation energy and vibrational frequency are compiled in Table 
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3. It is seen that excellent agreement with experiment (to 0.02 A, 0.1 eV, 
10 cm- 1 (l%)) is obtained at the (single-reference) CCSD(T) level. At the 
multi-reference level (without triples), the agreement is slightly less good, 
but this could certainly be improved upon by enlarging the 2-configuration 
reference space which was chosen in our calculations. 

3 Local increments 
3.1 Methodological aspects 

The main idea to be discussed here is the possibility to extract informa- 
tion from calculations on finite clusters and to transfer it to the infinite 
crystal. Such a transfer would certainly not be a good idea for global clus- 
ter properties (cohesive energy, ionization potential, etc.) - it only makes 
sense for local quantities. Now, localized orbitals are entities which can be 
defined within ionic crystals as well as within clusters of these materials. 
Moreover, electron correlation in or between such orbitals is a local effect. 
Therefore, if we prepare localized orbitals in the interior of a cluster (in a 
sufficiently solid-like environment) and if we calculate correlation energies 
involving these orbitals, we can hope to obtain transferable quantities. 

Of course, there is no reason to expect that these quantities would be 
additive in the solid. If we separately calculate, e.g., the pieces of correlation 
energy due to orbitals localized at ionic positions A, B, C, ... 



the correlation energy of the common orbital system of AB (or AC, BC, 
...) will deviate in general from the sum of constituents, due to inter- ionic 
interactions, and we can define non-additivity corrections such as 



Again, the next larger systems of three ions, ABC, will have correlation 
energies slightly different from the sum of the constituents plus the two-body 
non-additivity corrections, and this gives rise to three-body increments 



e(A),e(B),e(C),... 



(1) 



Ae(AB) = e(AB) - e{A) - e{B). 



(2) 



Ae(ABC) = e(ABC) - [e(A) + e{B) + e(C)] - 
[Ae(AB) + Ae(AC) + Ae(BC)}. 



(3) 



Similar definitions apply, in principle, to higher-body increments. 
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If we now make use of all these quantities, i.e., the intra-ionic correlation 
energies and the various inter-ionic correction terms, and multiply them with 
weight factors appropriate for the solid, we can hope to get a meaningful 
incremental expansion of the correlation energy per unit cell of the infinite 
crystal: 

6buik = E < A ) + ^ E A < AB ) + h E MABC) + ... (4) 

A A,B A,B,C 

In Ref. [4d] we have shown that this equation can be formally derived, under 
appropriate approximations, from an expression for the correlation energy 
of an infinite system. 

Let us discuss now the assumptions implicit in this approach more closely 
for the case of the MgO crystal. MgO is generally considered as a nearly 
perfect ionic crystal consisting of Mg 2+ and 2 ~ ions [17]; the question 
of a quantitative measure for the ionicity of MgO has been addressed only 
recently by Bagus and co-workers [17c, d], and in careful studies using various 
criteria the ionic charges have been shown to deviate from ±2 by <0.1 only. 
Thus, the attribution of localized orbitals to ionic positions made above 
seems to be a valid assumption. But even if there were some degree of 
covalency and/or some tendency for derealization in MgO, this would not 
invalidate our approach. In fact, the first applications of the method of 
local increments were made for covalently bonded crystals (diamond, silicon 
[4a,b]), and even for the 7r-system of graphite which according to usual 
classifications is considered as completely delocalized the method has been 
shown to yield meaningful results [4c]. 

Secondly, the determination of increments for non-additive inter-ionic 
correlation contributions in Eqs. 2 and 3 makes sense only, if the number 
of non-negligible increments is small, i.e., if the Ae(AB) rapidly decrease 
with increasing distance of the ions and if the three-body terms Ae(ABC) 
are significantly smaller than the two-body ones making the use of four- 
body contributions obsolete. A necessary pre-requisite for satisfying these 
conditions is the use of a size-extensive correlation method for calculating 
the increments. This excludes, for instance, the (variational) configuration- 
interaction method with single and double excitations (CISD), since it does 
not scale linearly with n for a system of non-interacting atoms A n . On 
the other hand, correlation methods of the coupled-cluster variety such as 
those discussed in the last section (ACPF, CCSD, CCSD(T)) do have this 
property. (These are benchmark methods of increasing complexity widely 
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used in quantum chemistry; we display results derived from all of them 
in the following Tables, in order to monitor convergence with respect to 
the many-particle basis set used.) When applying such a method, Ae(AB) 
should indeed rapidly decrease with increasing AB distance, since for non- 
overlapping pairs of ions only van der Waals-like correlation effects (~ 1/r 6 ) 
become effective. Furthermore, three-body terms can be expected to be 
significantly smaller than two-body terms, since two-electron excitations, 
involving a pair of orbitals at most, are known to dominate correlation 
effects. 

A final assumption underlying our approach is that of the transferability 
of localized orbitals from clusters to the bulk which was mentioned right at 
the beginning of this subsection. Its fulfillment depends, of course, on the 
preparation of the clusters. A free O 2- ion, e.g., would be unstable, and 
it is essential, therefore, to put this ion in a cage of Mg 2+ ions in order to 
stabilize it, and to simulate the Madelung potential of the surrounding ions 
in order to provide the correct field near the O nucleus. More details on 
cluster preparation and transferability tests will be given in the following 
subsections, where the determination of individual increments is discussed. 

3.2 Intra-ionic correlation 

The first increment to be calculated is the correlation energy which can be 
locally attributed to an 2 ~ ion in crystal surroundings of Mg 2+ and other 
2 ~ ions. As already mentioned, a realistic modelling of the crystal sur- 
roundings is essential, since otherwise the O 2- ion would not be stable at 
all. Fortunately, stabilization can be achieved in a both very simple and effi- 
cient way, by simulating the Pauli repulsion of the 6 nearest-neighbour Mg 2+ 
ions by means of pseudopotentials; we used the same energy-consistent pseu- 
dopotentials here as were used for the treatment of the Mg atom in the calcu- 
lations of Sect. 2. For representing the crystal environment of the resulting 
7-atom cluster, a Madelung approximation was made: 336 ions surrounding 
this cluster in a cube of 7x7x7 ions were simulated by point charges ±2 
(with charges at the surface planes/edges/corners reduced by factors 2/4/8, 
respectively). Here and in the following, the experimental bulk lattice con- 
stant (rMgO = 2.105 A) was adopted. Employing the [5s4p3d2f] basis set, 
which was already used for O and O - , for the description of the 2 ~ orbitals, 
too, we obtain the differential correlation-energy contributions, Ae(0 2_ ) = 
e(0 2 ~)-e(0), to the affinity of the extra electrons in crystal 2 ~ which are 
listed in Table 4. 
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A first point to make is that at all levels of approximation Ae(0 2 ) comes 
out considerably smaller than one would expect from a simple linear scaling 
of Ae(O-) values (2.77 vs. 1.86 eV, at the CCSD(T) level, cf. Table 1); such 
a linear scaling, which approximately works for the iso-electronic systems, 
Ae(Ne) and Ae(Ne~), [18] probably fails for O 2- due to the increased spac- 
ing of excited-state levels, when compressing the O 2- charge density in the 
(Mg 2+ )6 cage. When comparing individual Ae(0 2_ ) values in Table 4, we 
observe that in our single-reference calculations (active space 2s - 2p), the 
effect of single and double excitations is quite similar for ACPF and CCSD, 
while the inclusion of triples in CCSD(T) yields an increase by another 5%. 
Thus, the effect of triples is of less (relative) importance than in free O - but 
is still non-negligible. As in the case of 0~, we checked that enlarging the ac- 
tive space (to 2s - 3p) in the ACPF calculations (MRACPF) leads to a result 
numerically nearly identical to CCSD(T). Moreover, we tested the influence 
of an increase of the basis set ([5s4p3d2f] — > [6s5p4d3f2g] , cf. Table 1); at 
the CCSD(T) level, the correlation-energy increment changes by -0.008 a.u. 
(-0.2 eV), in line with our estimate given in Sect. 2. We also tried increasing 
the [5s4p3d2f] basis set by adding off-center functions (the (4s4p) sets of 
Sect. 2 at the positions of the Mg 2+ ions); the correlation-energy change 
of -0.005 a.u. is somewhat smaller here because not all components of the 
higher polarization functions (/, g) at the oxygen site cannot be simulated 
this way. A last test concerns the influence of the Madelung field. Leaving 
out all of the point-charges and performing the calculation with the bare 
(0 2 ~)(Mg 2+ )6 cluster leads to quite negligible correlation-energy changes of 
< 4 • 10 -5 a.u. only, at all theoretical levels; this underlines the notion of 
electron correlation being a local effect. Summarizing, lack of completeness 
of the one-particle basis set seems to be the largest source of error in the 
2 ~ results listed in Table 4, and the order of magnitude of the resulting 
error for the MgO cohesive energy per unit cell (with respect to separated 
neutral atoms) can be assessed to 0.2 ... 0.3 eV. 

Let us next consider the Mg — > Mg 2+ correlation-energy increment. This 
increment can directly be evaluated using the atomic calculations described 
in Sect. 2. This is so, since the Madelung effect identically vanishes here, 
when Mg 2+ is described by a bare pseudopotential. We checked this point 
by performing all-electron calculations for Mg 2+ with and without Madelung 
field: the non-frozen-core effect obtained thereby is of the order of 1 • 10~ 6 
a.u. with our largest basis set. Thus, the pseudopotential approximation is 
certainly valid here. 

Adding the Mg 2+ and O 2- correlation-energy increments of Table 4, we 
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get, at the highest theoretical level (CCSD(T)), a correlation contribution 
of 0.0547 a.u. to the bulk cohesive energy, £7 co h, per (primitive) unit cell 
of the MgO crystal. The experimental cohesive energy corrected for zero- 
point energy is known to be 0.3841 a.u. [19]; the most recent (and probably 
best) HF value is 0.2762 a.u. [6c]. This yields an 'experimental' correlation 
contribution to E co ^ of 0.108 a.u. which is just about the double of the intra- 
ionic value calculated so far. Thus, it is clear that inter- ionic contributions 
to be dealt with in the next subsection play an important role. 

3.3 Two-Body Corrections 

In this subsection, non-additivity corrections are determined, which arise 
when simultaneously correlating two ions in a cluster. 

Let us first consider here the interaction of ions with charges of opposite 
sign, i.e. Mg 2+ and 2 ~, which are next neighbours in the crystal. Using 
the same cluster as in the previous subsection when determining the intra- 
ionic correlation energy of 2 ~ (0 2 ~ plus 6 surrounding Mg 2+ plus 336 
point charges) and adding a core-polarization potential at one of the Mg 2+ 
neighbours of the central O 2- ion, we obtain an inter-ionic core-valence cor- 
relation contribution (cf. Table 5) which is due to the dynamic polarization 
of the Mg 2+ core by the O 2- valence electrons and which was clearly absent 
in the free Mg 2+ ion. Although the resulting value for the increment turns 
out to be considerably smaller than the intra-ionic correlation contributions 
of Table 4, its effect on the cohesive energy of MgO is by no means negligible, 
due to the large weight factor. 

In order to study the convergence of the correlation-energy increments 
with increasing distance of the ions, we replaced, in the cluster described 
above, one of the Madelung charges (at positions of Mg 2+ ions successively 
more distant from the central 2 ~ ion) by a pseudopotential and evaluated 
the influence of core-polarization. A rapid decrease with rM g o is observed, 
with the fourth-nearest neighbour MgO increment already approaching the 
numerical noise in our calculations. 

Turning now to the increments related to pairs of ions of the same kind, 
we can safely neglect Mg 2+ -Mg 2+ interactions. The correlation-energy con- 
tributions are exactly zero, at the pseudopotential plus core-polarization 
level. All-electron test calculations yield very small values around ~ 2 • 10~ 5 
a.u. for a nearest-neighbour pair of Mg 2+ ions in the crystal. 

Far more important are interactions between O 2- ions with their diffuse, 
fluctuating charge distributions. For the increment between nearest neigh- 
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bouring 2 ~ ions, we used a cluster with 448 ions, where two central O 2- 
ions were treated explicitly, while the 10 next Mg 2+ neighbours of these ions 
were simulated by pseudopotentials and the remaining ions were represented 
by point charges. The non-additivity correction with respect to two single 
2 ~ ions (cf. Table 5) turns out to be of the same order of magnitude as 
the MgO increment. The greater strength of interaction in the 0-0 case 
compared to the Mg-0 one (larger polarizability of 2 ~ compared to Mg 2+ ) 
is effectively compensated by the enhanced ion distance. 

Results for 2- -0 2 ~ increments between 2 nd through 5 th neighbours 
are also given in Table 5. The decrease with increasing ion distance is 
not quite as rapid as that for the Mg-0 increments. By multiplying the 
0-0 increments with r~ 6 , one can easily check that a van der Waals-like 
behaviour is approached for large r, and one can use the resulting van der 
Waals constant to obtain an estimate of the neglected increments beyond 
5 th neighbours; this estimate which is ~ 3 • 1CT 4 a.u. (including appropriate 
weight factors) should be considered as an error bar for the truncation of 
the incremental expansion of the MgO cohesive energy in our calculations. 

Adding up all the two-body increments which have been determined in 
this work, we find (cf. last row in Table 5) that the inter-ionic two-body 
correlation contribution to the cohesive energy of MgO is of comparable 
magnitude as the intra-ionic one. Thus, the important conclusion may be 
drawn that even in a (nearly purely) ionic crystal a mean-field description 
of inter-ionic interactions (electrostatic attraction, closed-shell repulsion) is 
not sufficient. 

3.4 Three-Body Corrections 

Let us complete now the information necessary for setting up the incremen- 
tal expansion of the correlation energy of bulk MgO, by evaluating the most 
important non-additivity corrections involving three ions. These corrections 
are obtained for triples with at least two pairs of ions being nearest neigh- 
bours of their respective species, using Eq. 3 (cf. Sect. 3.1). The numerical 
results are listed in Table 6. It is seen that the largest 3-body corrections 
are smaller by nearly two orders of magnitude compared to the leading two- 
body ones, thus indicating a rapid convergency of the many-body expansion 
with respect to the number of atoms included; the total three-body contri- 
bution to the correlation piece of the bulk cohesive energy, E co h is ~ 2% of 
the two-body part and of opposite sign. 
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3.5 Incremental Expansion 

In Table 7, the sum of local correlation-energy increments to the cohesive- 
energy, -Ecoh, of MgO (with respect to separated neutral atoms) is compared 
to the difference of experimental and HF values for E co ^. Our calculated 
values amount to between ~ 80 and ~ 85 %, depending on the correlation 
method applied, of the experimental value. The inclusion of triple excita- 
tions in the correlation method seems to be significant for describing the 
large fluctuating 2 ~ ions. A major part of the remaining discrepancy to 
experiment is probably due to deficiencies of the one-particle basis set: as 
discussed in Sect. 3.2, extension of the O 2- basis set to include g functions, 
for the evaluation of the intra-ionic contribution, already reduces the error 
by a factor of 2. 

A comparison to related theoretical results is possible at the density- 
functional level. A correlation-energy functional including gradient correc- 
tions (Perdew 91) yields a AE co h value of 0.087 a.u. [6b], only slightly inferior 
to our CCSD(T) one. Further density-functional results for lattice energies 
have been reported in Ref. [7]; these results have been determined for fixed 
HF densities, and the reference of the lattice energies is to Mg 2+ + O" + 
e~. As already mentioned in Sect. 1, the LDA value is much too high (by 
0.069 a.u.); the results with the Perdew-91 correlation-energy functional on 
top of the HF exchange is much better (error 0.014 a.u.), and the same may 
be said for a xc-DFT treatment with the Becke and Perdew-91 gradient cor- 
rections for exchange and correlation, respectively (error 0.011 a.u.). Using 
the Mg 2+ and O - data of Tables 1 and 2, our present calculations lead, at 
the CCSD(T) level, to a value of 0.071 a.u. for the correlation contribution 
to the lattice energy; adding this value to the HF result of Ref. [7], a de- 
viation from experiment of 0.009 a.u. is obtained, which is similar to that 
of the gradient-corrected DFT ones. Note, however, that this comparison 
must be taken with a grain of salt: the HF calculation of Ref. [7] is certainly 
not too accurate, since a relatively small basis set was used. Otherwise it 
would be difficult to explain why our correlation contribution to the lattice 
energy should be too large, although all possible sources of errors (discussed 
in Sects. 3.2 and 3.3) point to the opposite direction. 

4 Conclusion 

The correlation energy of the MgO crystal can be cast into a rapidly con- 
vergent expansion in terms of local increments which may be derived from 
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finite-cluster calculations. One-particle basis sets of triple- zeta quality with 
up to / functions at the positions of the O atoms and high-level quantum- 
chemical correlation methods (CCSD(T)) are necessary for obtaining ~ 85% 
of the correlation contribution to the bulk cohesive energy. About one half 
of this contribution can be attributed to intra-ionic interactions, the rest is 
due - to about equal parts - to dynamic polarization of the Mg 2+ cores by 
the 2 ~ ions and to van der Waals-like 2 ~-0 2 ~ interactions. 
Although numerical results from quantum-chemical calculations of this type 
are not necessarily superior to DFT ones, they provide additional physical 
insight into the sources of correlation contributions to solid-state properties. 
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Table 1: Electron affinity (eV) for the O atom at various theoretical lev- 
els (RHF: restricted Hartree-Fock self-consistent-field, CASSCF: complete 
active space self-consistent-field, (MR) CI: (multi-reference) configuration in- 
teraction, (MR)ACPF: (multi-reference) averaged coupled pair functional, 
CCSD(T) coupled-cluster with single and double substitutions (and per- 
turbative correction for triples); the reference configurations comprise all 
possible distributions of electrons within the space of active orbitals. 



basis set 


[5s4p3d2/H 


[6s5p4d3f2g] a 


[6s5p4d3/2 5 ] a 


active space 


2s - 2p 


2s - 2p 


2s-3p 


RHF 


-0.53 


-0.54 


-0.54 


CASSCF 


-0.53 


-0.53 


0.77 


(MR) CI 


1.02 


1.07 


1.36 


(MR)ACPF 


1.23 


1.28 


1.40 


CCSD 


1.18 


1.24 




CCSD(T) 


1.33 


1.40 




exptl. 


1.46 



"Correlation-consistent augmented polarized valence triple- and quadruple-zeta basis 
sets from Refs. [10,8a], contracted from (Ils6p3d2/) and (13s7p4d3/2i?) Gaussian primi- 
tives, respectively. 



Table 2: Ionization potentials Mg— >Mg + /Mg + — >Mg 2+ of the magnesium 
atom (in eV), from calculations using a 2- valence-electron pseudopotential 
(PP); in the (single-reference) ACPF calculations a core-polarization poten- 
tial accounts for core-valence correlation. 



basis set^ 


(4s4p) (4s4pld) 


RHF[ 

ACPF 


1 


6.63/14.75 6.63/14.75 
7.63/15.02 7.65/15.02 


exptl. 


7.65/15.04 



"Valence basis set optimized for the PP; exponents of Gaussian primitives are 
s: 2.4257193, 0.8226250, 0.1077490, 0.0394850; p: 0.7690470, 0.1886750, 0.0751010, 
0.0294970; d: 0.2; the basis set is used without contraction. 

6 For explanation of acronyms, cf. Table EL 
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Table 3: Bond length R e (A), dissociation energy D e (eV) and vibrational 
frequency u> e (cm -1 ) of the MgO molecule, evaluated at various theoretical 
levels. A Mg 2+ pseudopotential was used, and, at the correlated levels, a 
core-polarization potential was added. In the multi-reference calculations, a 
5-orbital (Mg 3s, O 2s, 2p) active space was chosen. 





Re D e U) e 


RHFQ 

MRCP 
MRACPF a 

CCSD a 
CCSD(T) a 


1.715 - 794 

1.750 2.01 751 

1.751 2.29 755 
1.744 1.96 746 
1.731 2.49 792 


exp-0 


1.749 2.61±0.22 785 





Tor explanation of the acronyms, cf. Table 1. Basis sets: O [5s4p3d2/]; Mg (4s4p). 
6 Ref. [16] 



Table 4: Intra-ionic correlation-energy increments and total contribution, 
AFj co h(MgO), of these increments to the cohesive energy of bulk MgO (in 
atomic units, 1 a.u. = 27.2114 eV). 





ACPF 


CCSD 


CCSD(T) 


Weight" 


Mg^Mg 2 +J 


0.04690 


0.04690 


0.04690 


1 


o^ o 2 - g 


-0.09548 


-0.09646 


-0.10162 


1 


AE coh (MgO)H 


0.04858 


0.04956 


0.05472 





"Weight factor in the incremental expansion of the bulk correlation energy (per primi- 
tive unit cell) of MgO. 

6 (4s4p)-valence basis set, cf. Table || PP+CPP description of Mg 2+ . 

c [5s4p3d2/] basis set, cf. Table |l[ Stabilization of O 2 by means of 6 Mg 2+ pseudopo- 
tentials at the nearest-neighbour positions of the MgO crystal (rM g o=2.105 A). Madelung 
field of outer ions in 7x7x7 cube represented by point charges ±2. 

d Correlation-energy change per primitive unit cell, for separation into neutral ground- 
state atoms. 
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Table 5: Inter-ionic two-body correlation-energy increments and total con- 
tribution, AE co j 1 (MgO), of these increments to the cohesive energy of bulk 
MgO (in a.u.). The notation A-B — > n means that the increment describes 
a nth nearest neighbour pair of A and B ions in the crystal. 



Mg-0 -> 1 ft 
Mg-0 -> 2 | 
Mg-0 -> 3 □ 




AE coh (MgO) 



ACPFflfl CCSD a ' b CCSD(T) a ' b 



-0.003031 
-0.000082 
-0.000014 



-0.003040 
-0.000083 
-0.000013 



-0.003037 
-0.000084 
-0.000014 



-0.002740 
-0.000238 
-0.000066 
-0.000027 
-0.000014 



-0.002538 
-0.000225 
-0.000061 
-0.000024 
-0.000012 



-0.002973 
-0.000264 
-0.000072 
-0.000029 
-0.000014 



0.037454 0.036139 



0.039066 



wci 



ghtf 



24 



6 
3 
12 
6 
12 



a For the definition of the increments, see Eq. 2. 

6 The (4s4p) valence basis set was used for Mg (Table ^) and the [5s4p3d2/] basis for 
O (Table]!]). 

c Weight factor in the incremental expansion of the bulk correlation energy (per primi- 
tive unit cell) of MgO. 

O 2- at (0,0,0) with 6 surrounding Mg 2+ carrying a PP (one of them an additional 
CPP) in a 7x7x7 cube of point charges. 

e Mg 2+ carrying a PP+CPP at (0,0,0), 2 ~ at (1,1,1) with a cage of 6 surrounding 
Mg 2+ PP, in a 7x7x7 cube of point charges. 

^same as e), but Mg 2+ at (-1,0,0) and O 2 " at (1,1,0). 

9 2_ ions at (0,0,0) and (0,1,1), 10 neighbouring Mg 2+ carrying a PP, remaining ions 
of a 7x8x8 cube represented by point charges. 

h 2 ~ ions at (±1,0,0), 11 Mg 2+ PP, in a cube of 9x9x9 ions. 
i 2 " ions at (0,0,0) and (2,1,1), 12 Mg 2+ PP, in a cube of 9x9x9 ions. 
j 2 ~ ions at (-1,-1,0) and (1,1,0), 12 Mg 2+ PP, in a cube of 9x9x9 ions. 
fe 2 " ions at (-2,0,0) and (1,1,0), 12 Mg 2+ PP, in a cube of 9x9x9 ions. 
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Table 6: Inter-ionic three-body correlation-energy increments and total con- 
tribution, AE co h(MgO), of these increments to the cohesive energy of bulk 
MgO (in a.u.). The notation A-B-C — * n means that AB and BC are 
nearest neighbours of their species, while AC are n th neighbours. 





ACPP^fQ CCSD a ' b CCSD(T) a ' 6 


weight^ 


0-0-0 -» 1 3 
0-0-0 -> 2 Q 


0.000069 0.000063 0.000067 
0.000014 0.000016 0.000013 


8 
12 


0-Mg-O -> 1 [J 
0-Mg-O ^20 


0.000013 0.000014 0.000016 
-0.000023 -0.000022 -0.000026 


12 
3 


AE coh (MgO) 


-0.000807 -0.000798 -0.000806 





a For the definition of the increments, see Eq. 3. 

6 The (4s4p) valence basis set was used for Mg (Table ^) and the [5s4p3d2/] basis for 
O (Table 0). 

c Weight factor in the incremental expansion of the bulk correlation energy (per primi- 
tive unit cell) of MgO. 

O 2- ions at positions (1,0,0), (0,1,0), (0,0,1), surrounded by a cage of Mg 2+ PP each, 
within a cube of 7x7x7 ions, with point charges ±2 simulating the outer ions. 

"Same as d), but 2 ~ positions (1,0,0), (-1,0,0), (0,0,1). 

/ Mg 2+ PP+CPP at (0,1,0), O 2- ions surrounded by Mg 2+ cage at (0,0,0) and (0,1,1) 
in a 7x8x8 cube of point charges. 

9 Mg 2+ PP+CPP at (0,0,0), O 2 " ions surrounded by Mg 2+ cage at (1,0,0) and (-1,0,0) 
in a 7x7x7 cube of point charges. 
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Table 7: Total correlation contribution, AE co h(MgO), to the bulk cohesive 
energy per primitive unit cell of MgO with respect to separated neutral 
atoms, from an expansion using the local increments of Tables 4 through 6, 
compared to the 'experimental' value. All data in a.u.. 





AE coh (MgO) 


ACPF 
CCSD 
CCSD(T) 


0.0852 (79%) 
0.0849 (79%) 
0.0930 (86%) 


exptl.0 
c-DFTQ 


0.1079 
0.087 



"Difference of the experimentai [19] and the HF value [6c] of the MgO cohesive energy; 
the experimental value has been corrected for zero-point energy. 

b Gradient-corrected correlation-energy density functional, Ref. [6b]. 
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